Quantum dynamics of two capacitively coupled superconducting islands via Josephson 

junctions 
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In this paper, we consider a system consisting of two capacitively coupled superconducting islands 
via Josephson junctions. We show that it can be reduced to two coupling harmonic oscillators under 
certain conditions, and solved exactly in terms of a displacing transformation, a beam-splitter-like 
transformation, and a squeezing transformation. It is found that the system evolves by a rotated- 
squeezed-coherent state when the system is initially in a coherent state. Quantum dynamics of the 
Cooper pairs in the two superconducting islands is investigated. It is shown that the number of the 
Cooper pairs in the two islands evolves periodically. 
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Recent experiments 1, 2] have demonstrated the pos- 
sibility of controUing or manipulating macroscopic quan- 
tum states in a single- Cooper-pair box, which consists of 
a small capacitance supuerconducting island coupled to a 
bulk superconductor via Josephson tunneling. Due to the 
small capacitance, there can be only zero or one excess 
Cooper pair on the island, which leads to an effective 
two-level system at appropriate values of external bias 
voltage. The possible use of these two-level systems as 
building blocks of quantum computer had been suggested 
already prior to the experimental work However, 
it is clear that realizing such a quantum computer even 
with a modest number of qubits will prove exceedingly 
difhcult and progress will be made only in small steps. A 
first step in this direction will consist in the coupling of 
two such Cooper-pair boxes. 

Marquardt and Bruder studied the quantum dy- 
namics of a system which consists of a Cooper-pair box 
and a capacitively coupling large superconducting island, 
predicted generation of mesoscopically distinct quantum 
states. The large superconducting island can be de- 
scribed by a harmonic oscillator since it has a compar- 
atively small charging energy, then we have Ej/Eq ap- 
proache the infinity, where Ej and Eq denote the Joseph- 
son coupling energy and the charging energy of the large 
island, respectively. The coupling system of a Cooper- 
pair box and a large superconducting island is analogous 
to a system which consists of a two-level atom interact- 
ing with a single mode of the quantized electromagnetical 
field. 

On the other hand, recently much attention has been 
paid to continuous-variable quantum information pro- 
cessing |llllll3|l3,|ll|lllilij which is based on har- 
monic oscillators in infinite dimensional Hilbert spaces. 
Therefore, it is is of significance to study such systems 



of capacitively coupling large superconducting islands 
via Josephson junctions in order to explore continuous- 
variable quantum information processing realizations in 
solid-state systems. In this paper, we consider theoret- 
ically another possibility: two large superconducting is- 
lands are capacitively coupled via two Josephson junc- 
tions. Each island may be described by a harmonic os- 
cillator, thus the system is reduced to a two-coupling- 
harmonic-oscillator system. We shall give an analytical 
solution of this system and investigate the quantum dy- 
namics of the Cooper pairs in the islands. 

Consider a system which consists of two capacitively 
coupled superconducting islands described in Fig. 1 
where Cj. and Cg. are the capacitances of z-th Josephson 
junction and i-th gate capacitor, respectively, Cc is the 
coupling capacitance, and Vi denotes gate voltage. The 
total energy of this system is the sum of the charging 
energy and the Josephson coupling energy. The charging 
energy is given by 
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where Qi is the total charge on the i-th island, Qc the 
charge on the coupling capacitor, and Qji a-nd Qg. de- 
note the charges on i-th Josephson junction and i-th gate 
capacitor, respectively. 

The charging energy can be expressed as a function of 
the total charges on the two superconducting islands Qi 
and Q2 alone and the circuit parameters. If let Ni — 
Qi/2e denote the number of the Cooper pairs on the i-th 
island with e being the charge of the electron, then the 
charging energy can written as 
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Ech = ^EcAN^-ng^f + Ei2{Ni-ng,){N2-ng^) 



'Ni — independent terms, 



(2) 



2 



where we have introduced the following notations 



Ei2 = 

Eci = 



CtiCt2 — 
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Ec2 — 



and the offsets introduced by the gates are given by 
ng^ = ^^[Cg^V.-{-irCc{V2-V,)\. 
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The Josephson coupling energy of the two superconduct- 
ing islands is given by 



Ejos = - ^ Ej^ cos ( 



(8) 
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Quantum mechanically, the number of the Cooper 
pairs and the Josephson phase are regarded as inter- 
conjugate operators, and satisfy the boson communica- 
tion relation [(j)j,Nk\ — iSjk- Then we get the Hamilto- 
nian 



H = 



Y.\EcAN^ 



Ej^ cos ( 



+Ei2{Ni-ng,){N2-ng2), 



(9) 



where we have discard the constant terms. 

For a large superconducting island, the Josephson cou- 
pling energy is much larger than the charging energy, i.e., 
EjjEci — + oo. In this case the approximation of a har- 
monic oscillator is valid, then one can replace cos (f>i term 
by the parabolic potential 1 — 0|/2. Hence the Hamilto- 
nian @ can be understood as that of two coupled har- 
monic oscillators. After making the following displacing 
transformation 



(10) 



the Hamiltonian @ can be expressed in terms of boson 
annihilation and creation operators as follows 
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Hi — hiOia\ai + \h{a\a2 + 02^1) 

1=1 

— A?i(a|a2 -f 0162), (11) 
where we have used the following decompositions 

1 _ (Ec^''' 
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(z = l,2), (12) 
(z = l,2), (13) 



FIG. 1: Circuit diagram for the considered system. The two 
superconducting islands appearing as boxes corresponding to 
the regions marked by dashed rectangles. The junction and 
gate capacitances are Cj. and Cg^ for i-island, respectively. 
The coupling capacitance is denoted by Cc and the gate volt- 
ages are Vi and V2. 



and the free-evolution frequencies and the coupling con- 
stant are given, respectively, by 



1 



= V^EciEj^ , 



(« = 1,2), 
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(14) 
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Then the displacement operator H1U|) can be rewritten 
in terms of the annihilation and creation operators as 



di{ng,,ng2) = Y[cxp 



i=l 



-i * - 



= £)(aoi,ao2), 
where agi and ao2 are given by 

1/4 



aoi = lUg^ 
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, ao2 = tng^ 



E, 



C2 



2E 



J2 
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(16) 
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For the sake of simplification, we consider the sym- 
metric situation where toi — uj2 = to since Ci — C2 and 
Cjj = . In this case, the cross term in the Hamil- 
tonian Hill) can be gotten rid of by the following beam- 
splitter-like transformation over 



B{(p) = exp (p{a[a2 — al^ai 



(18) 
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From and (|18|l we obtain the transformed Hamil- 
tonian 



+ f{ar+aj)-f{ol'+al), (19) 



where we have let lu[ = lui — X and Lu'2 = lu2 + A. In 
the derivation of the Hamiltonian (|19|) we have used the 
following formula 

W {Lp)aiB{ip) = di cos(/7 + 0,2 sin(p, (20) 
{ip)a2B{(p) — 02 cos</J — ai sincp. (21) 

Finally, we make a squeezing transformation over the 
Hamiltonian H19|l to get that 

i?3-^Uei,6)i?25(ei,6) (22) 
where the squeezing transformation is defined by 

^(6,6) = ^1(6)^2(6), (23) 

with Si{£,i) being the single- mode squeezing operator de- 
fined by 



'S'i(Ci) = exp 



e 



(j = l,2). (24) 



It is easy to show that when the squeezing parameters 
are chosen as 



(i = l,2), (25) 



the squeezing transformation can diagonalize the Hamil- 
tonian H2 as the following form 



^3 ~ hfliO^ai + 11^120202, 

where the frequencies il^ is given by 



(26) 



n,. 



A2 = 2e, 



El 



Cu-{~iycc' 



(* = 1,2) 



(27) 

In the derivation of the Hamiltonian (|23|l we have used 
the following formula 

Sl{£_^)a,Si{£,i) = flicosh^i - ajsinh^i, (28) 
sK^^)alS^{^^) = ajcosh^, -a^sinh^, (29) 

Assume that the two islands are initially in a coherent 
state 



|vI/(0)) =I)(ai,a2)|0,0), 



(30) 



where D[ai,a2) is the displacement operator with re- 
spect to the two modes defined by Eq. (|16|) . 



After making three anti-transformations <S''''(5i, ^2), 
B''{tt/A), and Z)'I'(q!oi, 002)7 the initial state becomes 

|*(0))3 = 5t(a,6)^^ (J) ^t(aoi,«o2)|ai,a2). (31) 

Then at an arbitrary time t, the wave function of the 
system is given by 



|*(i))3-exp(--i?3i) |*(0))3. 



(32) 



Making use of the following formula 

B^(J)i3(ai-aoi,a2-ao2)sQ) =^(71,72), (33) 

where the two parameters 71 and 72 are defined as 

7* - [(ai + (-l)^a2) - (aoi + (-1)^^02)] ,(34) 

We can express the evaluation of the state as 
2 

l*W>3 = n [^1 (-*M^(-e.)^(7j)] |0,0). (35) 

which indicates that when the system is initially in a co- 
herent state, the system evolves by a rotated-squeezed- 
coherent state in the transformed representation with a 
displacing transformation, a beam-splitter-like transfor- 
mation, and a squeezing transformation. It should be 
mentioned that the state (|35|l is an entangled state in 
the original representation, although it is a non-entangled 
state in the transformed representation. 

From (|32|l we can obtain the number of Cooper pairs 
on the i-th island 



Ni = "-gi + - [p2 sin(f]2i + S2) 



-{-iypism{n,t + Si)], (z = l,2). (36) 
with the amplitudes and phase shifts given by 



Pi = 



yiiEj 

tan"^ ( 



[{l^lif 



(Re7i 



,1/2 



~2C> 



tan ip.j^ ) 



(37) 
(38) 



Form Eq. H36|) it can be seen that the number of 
Cooper pairs on each island oscillates periodically, and 
one can control periods and amplitudes of these oscilla- 
tions through changing circuit parameters and gate volt- 
ages. 

The sum of Cooper pairs in the two islands is given by 

(A^i -I- iV2> - {ug, + ngj + p2 sin(f72t + 82), (39) 

and the difference of Cooper pairs in the two islands is 
given by 



(TVi - iV2) = (n. 



pi sin(f7i< + Si 



(40) 
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Making use of the wavefunction H35(l . we find that 
quantum fluctuation in each island is given by 

((AiV,)2) = (^gg-) ' ^[cos^ n,t + e'i^ sin^ J7,t],(41) 

which imphes that the quantum fluctuation of the Cooper 
pairs in each island oscillates periodically with the time 
evolution. The amplitude of the quantum fluctuation is 
dependent of the squeezing parameter ^i, i.e., the Joseph- 
son coupling energy and charging energy of the system 
under our consideration. 

In summary, we have investigated quantum dynamics 
of two capacitively coupled superconducting islands via 
Josephson junctions. We have shown that the system 



under our consideration can be reduced to two coupling 
harmonic oscillators. We have obtained an exact solution 
of the system when two harmonic oscillators are initially 
in coherent states, and found that the system evolves by 
a rotated-squeezed-coherent state. We have also stud- 
ied the Cooper-pair evolution and quantum fluctuations 
in the two superconducting islands and found that the 
number of the Cooper pairs in the two islands evolves 
periodically. 
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